Experimental study of the Ca 2 1 S+ 1 S asymptote 
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The filtered laser excitation technique was applied for measuring transition frequencies of the Ca2 
B-X system from asymptotic levels of the X 1 ^ ground state reaching v" = 38. That level has 
an outer classical turning point of about 20 A which is only 0.2 cm -1 below the molecular 1 S+ 1 S 
asymptote. Extensive analysis of the spectroscopic data, involving Monte Carlo simulation, allowed 
for a purely experimental determination of the long range parameters of the potential energy curve. 
The possible values of the s-wave scattering length could be limited to be between 250ao and lOOOao . 

PACS numbers: 31.50.Bc, 33.20.Kf, 33.20.Vq, 33.50.Dq 



I. INTRODUCTION 

In our previous paper [0 we reported on an accu- 
rate determination of the C&2 X x £+ ground state po- 
tential energy curve (PEC) from Laser Induced Fluo- 
rescense (LIF) spectroscopy. The study was motivated 
by the rapid progress in the development of the Ca 
frequency standard and the recent photoassociation 
spectroscopy data on the 1 So+ 1 Pi asymptote ||. It 
was considered as the first step of a comprehensive in- 
vestigation, which should provide a description of the X 
state for short and intermediate internuclear distances, 
thus enabling the precise determination of the long range 
part of the PEC, for example from photoassociation spec- 
troscopy or from a molecular beam experiment. 

The results of the LIF experiments, however, suggested 
that in the case of Ca2 it might be possible to make a 
reasonable description of the long range part of PEC al- 
ready by means of a classical spectroscopy. Here we mean 
the shallow potential well of the Ca dimer, which lead to 
thermal population even of the highest rovibrational lev- 
els for normal working temperatures (1220 K), so that a 
single laser scheme could be applied for their excitation. 
Of great importance is also the possibility to extend the 
description of the PEC by a long range dispersive expan- 
sion already from 9.4 A Q. 

For detecting transitions from ground state levels with 
v" > 35 we applied the Filtered Laser Excitation (FLE) 
technique, which was used also by Hofmann et al. B for 
studying the A 1 E+ excited state in Ca2- With this tech- 
nique we collected spectroscopic information for about 
24 highly excited levels of the X state with v" up to 38. 
The outer classical turning point of the last observed level 
with v" = 38 and J" — 10 is at 20 A and is located on 
the rotationally reduced potential only 0.2 cm -1 below 
the molecular asymptote. 

Although the FLE resolution was limited by the 
Doppler broadening, we will demonstrate that combining 
the data from [fij and this experiment, it is possible to 



reduce the uncertainty of the experimentally determined 
long range coefficients, and especially Cq to an extent, to 
compare them with the results of the most recent theo- 
retical predictions || |(| . 

In section [fl] we explain the experimental set up and 
discuss the achiev ed experimental resolution and uncer- 
tainties. Section III summarizes the methods for con- 



structing of the PEC. Analysis of the experimental data 
are performed in section IV , where we employed a Monte 
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Carlo simulation in order to determine the error limits of 
the fitted long range parameters of the PEC. 



II. EXPERIMENT 

The calcium dimers are obtained in a heat pipe oven 
already described in Ref. 0. The experiments are per- 
formed at helium or argon buffer gas under a typical 
pressure of 20 — 50 mbar and an oven temperature of 
1220 K. Working with argon as buffer gas enables us to 
reduce the pressure to a lower value than with helium 
without reducing the life time of the heat pipe due to 
solid calcium closing the optical path. In the frequency 
range accessible by the lasers in our group the most fa- 
vorable Franck-Condon factors between the levels close 
to the asymptote of the X X E+ state and the levels of 
the B 1 !]^" state (v' = 2) are less than 10~ 2 . In order 
to achieve a sufficient signal-to-noise ratio using Fourier 
transform spectroscopy a long term stability in the oper- 
ation of the oven is needed which is not achievable with 
our present apparatus. In addition, as it was mentioned 
in Ref.Jl]], the discrete fluorescence is accompanied by 
a strong background emission due to bound-bound-free 
and frec-bound-free molecular transitions. To overcome 
this problems we chose the filtered laser excitation (FLE) 
technique for direct excitation of the levels close to the 
asymptote of the ground state. 

Aim monochromator (GCA/McPherson Instru- 
ments) used as a narrow band pass filter of a typical 
window width 2 cm -1 is set on a strong line of the flu- 
orescence progression when exciting a transition to a se- 
lected rotational level with v' = 2 of the B state using a 
single mode dye laser (Coumarine 6). Then the laser op- 
erating around 550 nm (60 mW on the Ca2 sample) scans 
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FIG. 1: Experimental setup 



the frequencies of transitions from near asymptote levels 
of the ground state to the selected v' — 2 level. Laser 
frequencies resonant with such transitions will result in 
fluorescence, which matches the transmission band of the 
filter and can be detected. The signal after the monochro- 
mator is recorded as a function of the laser wavelength by 
a photomultiplier (Hamamatsu, R928). Although many 
other transitions of the X-B system contribute to the ab- 
sorption spectrum of C&2 in the scanned spectral region, 
only excitations which decay into the selected frequency 
window will be registered. Thus this selective technique 
provides greatly simplified absorption spectra in a region 
where the weak transitions from the asymptotic levels of 
interest are practically completely overlapped by other 
much stronger transitions. 

The relative positions of the FLE lines are calibrated 
using the frequency comb of a temperature stabilized 
Fabry-Perot interferometer with a free spectral range of 
149.75(1) MHz. The absolute line position is obtained 
using the absorption spectrum of iodine vapor at room 
temperature (fig. |l|) calibrated with the IodineSpec cal- 
culating software with which the positions of the iodine 
lines are predicted with an accuracy better than 3MHz 
in the studied frequency range 0, || . In order to improve 
the signal-to-noise ratio each spectrum is recorded sev- 
eral times. Before the averaging procedure, the recorded 
signal is normalized to the laser intensity to minimize the 
influence of its variations with frequency and time. The 
typical width of the observed lines is 0.05 cm -1 corre- 
sponding to the Doppler broadening. The intensity of the 
slowly varying background due to free-bound-free or free- 
bound-bound transitions of the B-X system is compara- 
ble to the intensity of the strongest lines and is mainly 
responsible for the noise in the spectra. Taking account 
of the linewidth, the signal-to-noise ratio and the uncer- 
tainty in the calibration with the I2 lines, the uncertainty 
of the absolute frequency for the strongest lines is esti- 
mated to be 0.0035 cm -1 . Two typical spectra are shown 
in fig]^ for low and intermediate J. The experimental 
traces as solid lines contain some gaps which were not 
scanned. The number above the line gives v" and the 
inset gives v' and J' of the excited state. 

Using the FLE technique we observed transitions from 
44 ground state levels and among them 24 close to the 
asymptote with v" > 35. The last observed levels are 
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FIG. 2: Portions of filtered laser excitation spectra obtained 
with different detection windows. As dotted lines the syn- 
thetic spectra used for assignment of the observed spectral 
features are shown. 



(v" = 38, J" = 6,8, 10) as summarized in the data field 
in fig. |3|. The assignment of the new observed lines is 
performed using the potential energy curve determined 
in our previous study of the ground state Q and a nu- 
merical potential of the B state determined from the data 
available in the literature |j] . Simulations of the transi- 
tions by these potentials for selected detection windows 
are presented in figure ^ as dotted lines. These lines are 
shifted vertically with respect to the observed traces in 
order to show clearly the quality of the simulation. The 
intensity of the synthetic spectrum is determined by the 
appropriate Franck-Condon factors, by the thermal pop- 
ulation of the ground state levels and by the selected 
detection window. The small differences between the ob- 
served and the simulated spectral features could be at- 
tributed to the limited accuracy of the B state PEC (0.05 
cm -1 ) and to the fact that the used ground state poten- 
tial was determined for levels with only v" < 35. Also, 
only transitions in the main calcium isotopomer, C&2 
were simulated. The unambiguity of the assignment was 
usually proved further by the observation of relatively 
long (5-6 vibrational quanta) and self consistent vibra- 
tional progressions of P,R doublets. 

With the present set-up we can study the influence of 
the collisions between calcium dimers and the buffer gas 
by changing its pressure since no additional instrumental 
broadening on the lines is present. We found that the 
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FIG. 3: The range of v" and J" quantum numbers of the 
observed levels in the X 1 Ej" state of 40 Ca2. 



III. CONSTRUCTION OF PEC 

Following the approach adopted in the previous report 
|jj , information which concerns only the ground state was 
extracted from the raw spectroscopic data. From transi- 
tion frequencies with a common upper state level, differ- 
ences between ground state levels were calculated. The 
PEC for the ground state was then constructed, which 
describes these experimentally observed differences. Sim- 
ilar to Ref.|IJ the possible combinations of line frequen- 
cies within a progression were restricted only to pairs of 
one P and one R lines, forming a total of 8500 differences. 

The methods for construction of PEC were discussed 
in detail in Ref. [Q and will be summarized only briefly 
here. 

The first method represents the potential as a trun- 
cated expansion over analytic functions: 



collisional broadening and the shift of the observed lines 
are much smaller than what was previously assumed [Q. 
This means that the experimental uncertainties of the 
transitions obtained in that study were overestimated. 

For completeness we enriched our data with progres- 
sions obtained by LIF using the 501 nm, 488 nm and 457 
nm lines of an Argon ion laser (Spectra Physics, Beam- 
Lok 2060) recorded by Fourier transform spectroscopy. 
Most of the transitions which are known to be excited 
by these lines [[l0| were observed in this study. The 
present set of transitions representing all the LIF and 
FLE data covers a range of rotational quantum numbers 
from J" — 4 to 164 and a range of vibrational quantum 
numbers from v" = to 38. We have recorded a total 
of 3580 transitions resulting in 924 levels of the ground 
state |§. 

During this study six atomic calcium lines were 
also observed around 1.9/im, which correspond to the 
3 Di,2,3 — > 3 Pi,2,3 transitions, when irradiating the sam- 
ple by Ar + laser lines. These transitions have been al- 
ready observed in Mg in a similar experiment JTT[ | . Here 
they could be attributed to predissociation of the molecu- 
lar B state through a coupling with repulsive states corre- 
lated to the 3 D + X S asymptotes. It is worth mentioning 
that the 3 Di2,3 — > 3 Pi,2,3 transitions were observed 
also at temperatures down to 950 K, which suggests that 
the first step of the process is predominantly excitation 
of pairs of free Ca atoms by the laser radiation to the 
molecular B state. The intensity ratio between the in- 
frared emission and the laser light was highest for the 
457 nm Ar + line and decreases with the increase of the 
laser wavelength. 

After a new examination of the previous Fourier trans- 
form spectra using the new potential energy curve which 
will be presented in the following, three unassigned weak 
molecular progressions were found to belong to the cal- 
cium isotopomer 44 Ca 40 Ca. The exciting transitions were 
identified as (0, 42)<-(4, 43) and (0, 72)^(5, 73) at 
18788.36 cm" 1 and (0, lll)<-(8, 112) at 18787.36 cm" 1 . 



where at, b and R m are parameters (R m is close to the 
equilibrium distance). For short and long intcrnuclcar 
distances the potential is smoothly extended respectively 
for R < i?i nn with 

A + B/R 12 (2) 

and for R > R ou t with 

D c - C 6 /R 6 - C s /R 8 - C w /R 10 . (3) 

Here D c is the value of the dissociation energy defined 
with respect to the minimum of the PEC. Since i? ut is 
usually chosen within the region where the contribution 
of the exchange energy is negligible, no additional dump- 
ing (or cut off) functions for the dispersion coefficients 
are introduced. 

Parameters and Cq and C$ are fitted in a nonlinear 
fitting procedure (for details see [fl2)), while A, B, C±o 
and D e are adjusted by the program in order to ensure 
smooth connection between the extensions and the an- 
alytic form. R m , b and the connecting points R uln and 
i? ou t are kept fixed to values, which allow fast conver- 
gence of the fitting routine. 

The second method defines the PEC for short and in- 
termediate internuclear distances as a set of points con- 
nected with cubic spline function. For large internuclear 
distances the same long range expansion as in (^) is used. 
The values of the potential Ui in a preselected grid of 
points, D c , Cq, Cg and Cio are treated as fitting param- 
eters. They are determined in an iterative fitting routine, 
which linearizes the problem by using a modified version 
of the inverted perturbation approach Jl3| . The connect- 
ing point i?out is chosen in the following way. Initially, 
the PEC is constructed in a pointwise form up to 13 A. 
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Then only the long range parameters are fitted assuming 
that i? ou t = 9.5 A. Although both sections of the PEC are 
determined independently, their shapes turn out to be al- 
most identical between 9.5 A and 13 A. As a last step the 
Cio coefficient of the long range expansion is slightly ad- 
justed (the change does not exceed few percent) in order 
to fit best the shape of the pointwise potential between 
9.4 A and 10 A. The crossing point of the pointwise and 
the long range curves is taken as i? ut- The new value of 
i? out differs from that which was used when fitting the 
long range parameters. Usually the change is small (typ- 
ically 0.1 A) and since the two sections are very close in 
this region, the effect of such a change on the quality of 
the fit is negligible. Of course, a second iteration could 
be performed with the new value of i? ut , but usually it 
is necessary only in the beginning, when the initial values 
of the fitting parameters are far from the best ones. 

Before fitting the PEC to the new data set, the exper- 
imental errors of the data obtained by LIF were reana- 
lyzed, since the value of the normalized standard devia- 
tion a ~ 0.45 given in the previous study (Ref. p]) signals 
their overestimation. The value of the uncertainty before 
was chosen to be 0.01 cm -1 for a strong LIF line in order 
to take into account possible frequency shifts due to col- 
lisions. Since the new FLE measurements showed that 
the probable influence of the temperature and the buffer 
gas pressure are much smaller than expected we reduced 
the errors of all LIF frequencies by a factor of 0.6, which 
will give a more realistic estimate on the error limits of 
the fitted PEC parameters. 

The analytic potential listed in table § describes the 
differences between the observed spectral lines from both 
experiments with a standard deviation a = 0.0064 cm -1 
and normalized standard deviation ct=0.69. The quality 
of the PEC for the near asymptotic levels is estimated 
by calculation the standard deviation with a reduced set 
of differences, where for each difference at least one level 
belongs to v" > 35. We obtained cr 35 = 0.00 92 cm" 1 
and = 0.92. The same parameters for the numeri- 
cal potential (table 0) are a = 0.0068 cm' 1 , a = 0.74 
and (T35 = 0.0092 cm -1 , a 35 = 0.89. In order to cal- 
culate the value of the pointwise potential in the range 
R < 9.44 A, a natural cubic spline through all the grid 
points should be used. The parameters of both curves 
are chosen such, that their minima are set to zero. In or- 
der to facilitate the comparison between the dissociation 
energies Dq defined with respect to the lowest rovibra- 
tional level (?/' = 0, J" = 0), their values are given in 
tables I and |l| as additional parameters. 



IV. UNCERTAINTIES OF THE LONG RANGE 
PARAMETERS 

The motivation of this study is not only to describe the 
rovibrational structure of the 40 Ca2 ground state. We 
rather want to examine to which extent the performed 
spectroscopic study is able to go beyond the reproduction 



TABLE I: Parameters of the analytic representation of the 
X 1 E^" state potential energy curve in 40 Ca2. 



R < 3.66 A 
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B 



3.66 A 

-2.9714 x 10 2 cm -1 
7.209 x 10 9 cm^A 12 
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R > 9.5 A 



Rout 

D c 
C 6 
C 8 
Cio 



9.5 A 
1102.076 cm -1 
1.0030 x 10 7 cm" 1 ! 6 
3.87 x 10 s cm _1 A 8 
4.39 x 10 9 cm" 1 A 10 



Additional parameter 
Do = 1069.871 cm" 1 



of the observed differences between ground state levels. 
Of greatest interest is the reliability of the determined 
long range parameters (and especially C 6 ) since they play 
an important role in phenomena involving two interacting 
cold Ca atoms in the ground state. 

Our analysis relies on two main assumptions: 

• the interaction between two Ca atoms can be de- 
scribed within an adiabatic picture, through a sin- 
gle channel model applying the one dimensional 
Schrodinger equation with an effective potential en- 
ergy curve; 

• the shape of the potential for internuclear distances 
R > 9.4 A can be described with the long range 
expression (E|), neglecting the exchange energy. 
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TABLE II: 
X X E+ state 



Parameters of the numeric representation of the 
potential energy curve in 40 Ca2- 
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5.809524 


684 


9589 


3 


280470 


4525 


7282 


5.940476 


728 


9235 


3 


372215 


3090 


9557 
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892857 
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8199 


9.076923 


1074 


5969 


5 


023809 


322 


3873 


9.435897 


1080 


8961 


5 


154762 


393 


7222 


9.794872 


1085 


5974 


5 


285714 


461 


4555 


10.303419 


1090 


2990 


5 


416667 


524 


6311 


10.811966 


1093 


5160 


5 


547619 


582 


9870 


11.611111 


1096 


6870 



D c = 1102.060 cm' 1 
Rout = 9.44 A 
C 6 =1.0023 xlO 7 cm -1 A 6 

Additional parameter 

D = 1069.868 cm" 1 



C 8 =3.808xl0 8 cm -1 A 8 
Ci =5.06-10 9 cm^A 10 



Due to the zero net electronic and nuclear spin in the 
case of 40 Ca and the large energy separation between 
the ground state and the lowest excited states the first 
assumption seems to be well justified. The second one is 
supported by the expected small value of the exchange 
energy jj], |lj] at 9.5 A and the dominant van der Waals 
character of the interaction. 

Under these assumptions we ask the question: what are 
the variations of the PEC parameters around the fitted 
ones which are still in agreement with the experimental 
data? 

In Ref. |jj the answer was given by plotting a contour 
plot of the likehood function % 2 by varying the long range 
parameters. For the nonlinear fitting procedures it is, 
however, not straightforward how to use the contour plot 
in order to determine the confidence interval of the pa- 
rameters since it is not obvious to which confidence limit 
a given contour corresponds. That is why there is some 
ambiguity in the probable errors of the fitted parameters 
determined only by using contour plots. The use of the 
matrix of variances and covariances, which is obtained 
by the fitting routines in the linearized form is not useful 



TABLE III: Parameters of the long range expansion for the 
X 1 Ej state in 40 Ca2 derived in this study and compared with 
the most recent data from the literature. 

This study Other sources 



D c , cm 1 

Do, cm -1 

Ce x 10 7 cm" 1 ! 6 



C 8 x 10 8 cm" 1 A 8 
C 10 x 10 9 cm" 1 ! 10 



1069.868(10) 
1.003(33) 



3.15 + 4.46 



1.7-7-8.4 



1095.0(5) (10 1 
1102.08(9) TI| 
1069.88(9) [1| 
1.070(6) [5| 
1.098 if 
1.15 no 
1.02 + 1.12 | 
3.27 |B 
1.1 + 3.8 |1 
4.74 |S 
3.7 + 17.0 j! 



in our case. The reason is that even if we assume the 
errors of the experimental frequencies being independent 
and following a normal distribution, the quantities which 
enter the fit are actually differences between these fre- 
quencies. Consequently their errors are not independent 
and we are not allowed to interpret the above mentioned 
matrix as in the case of independent data points. 

A possible way to get the probability distribution of the 
fitted parameters is to perform a Monte Carlo simulation 
of synthetic data. 



A. Monte Carlo simulation 

Since our experimental data come with some measure- 
ment uncertainties, which we assume to be normally dis- 
tributed, the set c(°) of the fitted parameters of the PEC 
will differ from the true one c^K If we perform a series 
of similar experiments the experimental data sets will be 
slightly different and consequently we will obtain differ- 
ent sets of fitted parameters cilp- Having a sufficient 
number of measurement sets we could plot the distribu- 
te) 

tion of c e xp and then decide what are the most probable 
values and what are their uncertainties. 

So let us assume that the fitted potential U^(R), 
described with a set of parameters c' ' , is not too far 
from the true one. Then, knowing the experimental un- 
certainty associated with each observed transition fre- 
quency, we can synthesize a set of experimental data 
which will have exactly the same structure as the real 
one (measurement errors, distribution of v and J quan- 
tum numbers). For this goal transition frequencies are 
formed using the eigenvalues for (R) (since the up- 
per state levels play no role in our analysis their ener- 
gies could be set to zero) and adding to each calculated 
frequency a small random quantity normally distributed 
with mean value zero and a standard deviation equal to 
the experimental error. Then we transform the synthetic 
frequencies into synthetic differences and use the same 
fitting procedure in order to obtain a new set of param- 



(i) 

eters c sim . Performing a large number of simulations we 



6 



can plot the distribution of c^ m — . Since the simu- 
lated data are related statistically to (R) in the same 
way as the real experimental data to the true potential 
and since these two potentials were assumed to be close 
to each other the obtained distribution should be not too 
different from that of ci% — . 

In order to realize these ideas we used the pointwisc 
representation of the potential since in the present real- 
ization of our computer codes it allows more flexibility in 
connecting the long range expansion with the potential 
curve at intermediate internuclear distances. The prize to 
pay for this is that the first derivative of the PEC might 
have a small "kink" in i? ut- This problem, however, 
does not influence the analysis of the PEC and it could 
be avoided by slight adjustment of the pointwise curve. 
What we win is that all the long range parameters are fit- 
ted to the experimental data independently of the actual 
shape of the potential at intermediate internuclear dis- 
tances. The continuity of the potential is ensured mainly 
by varying the connecting point i? ut and to a smaller ex- 
tend Ciq. Analysis similar to those, performed in Ref. jlj 
showed that the present body of experimental data fixes 
the inner part of the PEC approximately between 3.5 A 
and 9.5 A. Possible variations in this region are connected 
mainly with the experimental uncertainties and the in- 
fluence of the selected representation for the potential is 
negligible. Therefore we may expect that the choice of 
the connecting point i? ut around 9.5 A makes our anal- 
ysis almost independent of the functions used to model 
the PEC at short internuclear distances. 

Having generated a set of synthetic differences 

A-E&i.t** = E vhi ~ E V2J2 ^ treat them aS eX P er " 

imental data. U^°\R) is taken as initial potential and a 
small correction to it is fitted which should minimize the 
difference between A2?„, r „ , and the differences cal- 
culated with the corrected potential. The correction is 
defined as follows: 



N 

6U(R) = J2 SuiSi(R), for R < R out (4) 

2=1 

and 



SU(R) = 6D e - ^ - ^ - for R > R out , (5) 

where Si(R) are defining functions of the pointwise po- 
tential (see |l5[), 5ui are the values of the correction in 
an equidistant grid of internuclear distances Ri and SD e , 
SCe, SCg, SCw are the corrections to the long range pa- 
rameters. Using the theory of perturbations the shift of 
the difference A.E^j v j 2 due to 6U(R) can be written 



as: 

»=i 

+ 5D e (L v ^ h - L v ^) 
- 5Cj(Lf Jl - Lf J *) 

.7=6,8,10 



where 

Kf= [ R " Ut *%j{R)Si{R)dR (7) 
Jo 



Lf= / V 2 vJ (R)R- j dR (8) 

J Rout 

are the corresponding mean values of Si(R) and R~° 
calculated with the wave functions ^> v j(R) of the levels 
forming the difference. 

Note that in principle the exact presentation of U 1 - ^ (R) 
is not of importance since we study only the possible 
small variations around it which are still allowed by the 
experimental data. The advantage of expressing the cor- 
rection to the potential in a pointwise form (cubic spline 
function is used for interpolation) is that the influence of 
the parameters Sui is localized in R, since away from Ri 
Si(R) decays exponentially Q. This gives a high flexi- 
bility in constructing the needed form of the correction. 
Additionally, the correlations between Sui are introduced 
mainly by the experimental data and not by the selected 
basis functions. 

With each set of simulated data a new fit is performed 
which adjusts 5ui, SD e , SCq, 5C$ and SCiq. Along with 
this we examine the contributions due to Sui and the long 
range corrections to the fitted differences (the sums in 
Eq. ^) . As it could be expected, the changes of the PEC 
for R < R ou t axe very small (8ui ~ 0.01 cm -1 ) for all sets 
of synthetic data due to the large amount of experimental 
observations in this region. In addition the corrections to 
the differences due to Sui, i.e. the first term of the right 
side of eq. [| turn out to be smaller than the correspond- 
ing experimental uncertainties and in principle one could 
neglect them. We checked this and, indeed, the distri- 
bution of the long range parameters was almost identical 
for simulations with and without fitting Sui. Since this 
is true for different numbers of Sui parameters, we may 
consider it as a proof of the statement made above, that 
in the present case the determination of the long range 
parameters is almost independent on the model used to 
describe the inner part of the PEC. 

The positions of the energy levels and the molecular 
asymptote D e in Ref. [Q were defined with respect to 
the minimum of the PEC, U m i n . Then we estimated that 
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FIG. 4: Projections of the long range parameters distribution realized with 3000 simulated points on the Do-Ce and Ce-C's 
planes. CL stays for confidence limit. 



the possible variation of t/ m i n for different representations 
of the potential is of the order of 0.01 cm -1 , which was 
smaller than the uncertainty of the dissociation energy, 
0.09 cm -1 . With the extended experimental data set 
and after the revision of the experimental uncertainties 
we are now able to give a much better prediction of the 
dissociation energy with uncertainty also of the order of 
0.01 cm^ 1 . Obviously denning the dissociation energy 
with respect to U m { n will introduce an undesired increase 
of the uncertainty for this long range parameter. Con- 
sequently, starting from here we will consider only the 
value of the dissociation energy Do denned with respect 
to the lowest rovibrational level Eqq — -B„" = o,j"=o- In- 
deed, fitting Sui to the simulated data we convinced our- 
selves that although the shape of the potential (and also 
Umin) could vary slightly from fit to fit, the overall effect 
on the position of Eqq is much smaller than the typical 
variation of C/ m i n . Nevertheless, the contribution of Eqq 
to the uncertainty of Dq is taken into account. 

In this way, the main steps of the performed Monte 
Carlo analysis are: 

1. Find the PEC C/ (0) (i?) which can reproduce the 
experimental observations with smallest standard 
deviation; 

2. From the eigenvalues for (R) generate a set of 
synthetic data by adding normally distributed ran- 
dom deviations according the experimental uncer- 
tainties. This set has exactly the same structure as 
the experimental data; 

3. Perform a fit with the synthetic data as it was done 
with the original one and collect the obtained po- 
tential parameters in a distribution list; 

4. Repeat steps (2)- (4) until enough simulated long 
range parameters are collected in order to analyze 
their distribution. 



B. Results 

In figure ^ two projections of the long range parameters 
distribution obtained from 3000 simulations are shown. 
It is worth mentioning the strong correlation between Cg 
and Cg, which is the case also for the other pairs of dis- 
persive coefficients. In order to determine the confidence 
regions for the parameters we construct an ellipse in the 
four dimensional space of the parameters, whose axes 
have such orientations and relative lengths that the re- 
sulting surface follows the form of the distribution. Then 
we change gradually the lengths of the axes by a constant 
factor and count the number of points which are enclosed 
by the ellipse. For example, the ellipse which contains 
68.3 % of the points will determine the 68.3 % (la) con- 
fidence region of the parameters. The projections of the 
68.3 % and the 90 % confidence regions (strictly speaking 
the contours surrounding the projections) are plotted in 
the left part of figure [| for the pair of parameters Dq and 
C 6 . 

This allows us to determine the mean value and un- 
certainty of the dissociation energy of the 40 Ca2 ground 
state with respect to the lowest rovibrational level as 
1069.868 ± 0.010 cm" 1 and of the C 6 coefficient as 
(1.003 ± 0.033) ■ 10 7 cm" 1 ! 6 at a confidence limit of 
68.3 %. The values and the uncertainties of all the long 
range parameters are summarized in table III and com- 
pared with the values from other sources. 

The highest observed rovibrational level in Ref. [Q had 
v" = 35 and it was concluded that with the available 
experimental data set it is not possible to give any rea- 
sonable estimation for C§. Now, by using a Monte Carlo 
analysis we can be more rigorous. In figure [s] we compare 
the projections of two distributions on the Dq-Cq plane. 
The first one is derived using only the data from Ref. Q| 
before the revision of those experimental errors. In order 
to obtain a meaningful estimation of the dissociation en- 
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FIG. 5: Projections of the long range parameters distributions 
on the Do-C(j plane derived from the LIF data from Ref. jj]] 
and from this study (the black spot in the centre). 



FIG. 6: Projections of the long range parameters distributions 
on the Dq-C§ plane derived by assuming the validity of the 
long range expansion from 9.4 A (solid circles) and from 10 
A (open circles). 



ergy of the X state, we restricted in Ref. jD] the possible 
variation of Gq within a ± 5 % interval around the the- 
oretical value ||. This interval is denoted in fig. || with 
horizontal dashed lines and leads to an uncertainty of the 
dissociation energy similar to that, given in Ref. §. Af- 
ter including the new LIF and FLE data the distribution 
shrinks to the small black spot in the centre of figure || 
(see also the inset), which is actually the distribution 
from figure |]. It would be, however, not correct to inter- 
pret the drastic reduction of the size of the distribution 
only as a result of the transitions involving rovibrational 
levels with v" > 35. In fact the items, which led to the 
distribution presented in figure ^ are: 

• The assumption that the long range expansion (|^) 
is valid starting from 9.4 A. This allows to link the 
values of the long range parameters with the fre- 
quency differences involving levels with 30 < v" < 
38, i.e. the shape of the long range potential is 
roughly tested from 9.4 A up to w 20 A (the clas- 
sical turning point of the last observed level). A 
Monte Carlo simulation shows (figure ^|) the change 
of the distribution when shifting i? out from 9.4 A 
only by 0.6 A to i? out = 10 A. Note, that in this 
case the loss of accuracy on C§ is larger by a factor 
of two compared to the increase of the Do uncer- 
tainty; 

• Availability of highly excited levels with v" = 38, 
since they limit the variation of the dissociation 
energy; 

• The enlargement of the experimental data, which 
allow accurate determination of the PEC for short 
and intermediate internuclear distances. Because of 
the large body of data, the long range parameters 
are almost independent of the exact representation 



of the rest part of the PEC, which is clearly indi- 
cated by the small variations of Sui (eq. ||) during 
the statistical analysis; 

• The corrected estimate of the experimental uncer- 
tainties of the data in Ref. [Q . 



C. The value of the Cs coefficient 

The confidence interval for Cq given in this paper 
(1.003 ± 0.033) • 10 7 cm^A 6 does not overlap with the 
most recent theoretical prediction of (1.070 ± 0.006) • 10 7 
cm _1 A 6 within the error stated in Ref. ||. 

From our side, we realize that an important point in 
our analysis is the assumption for the validity of the long 
range expansion (§) down to w 9.4 A. As it was shown in 
fig. shifting R out to larger internuclear distances will 
increase the uncertainty on the fitted parameters. The 
same will happen if we introduce additional terms into 
the long range model (for example the exchange energy) . 
Our analysis can readily include these additional param- 
eters, provided there are new experimental observations 
which require this. 

D. The s-wave scattering length 

In order to obtain the s-wave scattering length and its 
confidence interval, it is straightforward to calculate its 
value for each of the simulated PEC. These calculations 
will show how reliably the spcctroscopically determined 
potential up to « 20 A can be applied to model cold 
collisions. The results are presented in fig. frla. The po- 
sition of the last bound level (v" = 40) with respect to 
the molecular assymptote, Dq — E v >r = 4o,j"=o is shown in 
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fig. 0b. The values of the scattering length a within the 
68.3 % confidence region vary from 250ao to f 000a (the 
Bohr radius a « 0.52918 A). Since these values indicate 
a bound level very close to the asymptote (several MHz, 
see fig. ^b) , initially we were surprised that from our much 
less accurate data (several 100 MHz) we can make such 
a good prediction of the position of the last bound level. 
Obviously the large amount of observed data fixes the ac- 
cumulated phase of the wave function for this level along 
the PEC up to w 20 A and the variation of the binding 
energy with respect to the asymptote is restricted mainly 
by the uncertainty of Cq. In addition, from figure 0b we 
see that the smaller the binding energy, the less sensitive 
it is to the uncertainties of the other long range param- 
eters. 

The inner wall of the potential (R < 3.5 A), however, is 
not well fixed by the experimental data. Could its change 
shift the last level by several MHz without influencing 
the positions of the other levels? The answer is negative. 
The wave functions for small intcrnuclcar distances for 
the highly excited levels differ only by a constant factor. 
For example, for a PEC with binding energy Dq — E40 « 
3 MHz for v" = 40, the shapes of the wave functions for 
v" = 40 and v" = 35 are almost identical up to 7.5 A 
and the ratio between them is: 



35 for R < 7.5 A. 



(9) 



So, using the perturbation theory it is clear, that a 
change of the potential for R < 3.5 A, which will shift 
the position of the level with v" = 40 by 1 MHz will 
shift the level with v" = 35 by « 1.2 GHz (0.04 cm" 1 ). 
Since this change will not alter for example the position 
of the lowest level (v" — 0), the introduced shift of the 
difference £v'=35 — E v n— will exceed the experimental 
uncertainty roughly by a factor of 6, which is a contra- 
diction. 

Due to the revised confidence interval for the Cg, the 
uncertainty of the scattering length was not reduced 
compared to Ref. [Q (112a -j- 800ao), although the un- 
certainty of Cq itself is smaller than the assumed one 
in Ref. JpJ (±5 % around the theoretical prediction of 
1.07 • 10 cm _1 A 6 ||). Following the dependence of the 
binding energy on the last level on Cg (fig. 0b) we see 
that for smaller values of Cq this level approaches the 
asymptote, which makes an accurate determination of 
the scattering length more difficult. 



V. CONCLUSIONS 

The 1 S+ 1 S molecular asymptote of the Ca dimer was 
studied by employing the filtered laser excitation tech- 
nique. The spectroscopic data obtained in our previous 



study by Fourier transform spectroscopy were enriched 
by adding 56 transition frequencies from ground state 
levels with v" up to 38. This gave us an opportunity to 
reanalyze the existing description of the long range part 
of the PEC from Ref. §. 

By assuming the long range expansion ([|) for the po- 
tential energy curve to be valid already from 9.4 A, the 
probability distributions of the dissociation energy and 
the dispersion coefficients were derived from Monte Carlo 
analysis. Their mean values and uncertainties are given 
in Table EH. We showed, that due to the large amount 



of experimental data, the long range analysis is almost 
independent of the model functions used to describe the 
inner part of the PEC. 

The reliability of the experimental potential for de- 
scribing collisions between two Ca atoms at low temper- 
atures was checked by calculating the s-wave scattering 
length and the position of the last bound level with re- 
spect to the asymptote for the variety of long range ex- 
tensions of the potential allowed by the present experi- 
mental data. Following the distribution of the long range 
parameters, the binding energy of the last level was found 
to vary between several 10 kHz and several MHz. This 
confirmed our anticipations 0] for a large and positive 
scattering length, 250ao < a < 1000ao. Contrary to 
the previous study, the present determination is purely 
experimental and does not rely on the theoretical esti- 
mations for Gq. Moreover, the difference in the values of 
C6 from Ref. S and from this study exceeds the stated 
uncertainties, although a small model dependence is still 
existing by altering the connecting point between the in- 
ner potential and the long range branch (compare fig. |^). 

The Monte Carlo simulations open new perspectives 
for analyzing the uncertainties of the parameters of the 
PEC. They could be used also to determine the kind of 
spectroscopic data and their required accuracy for achiev- 
ing the desired uncertainty of the fitted parameters. For 
example, we checked that including a vibrational pro- 
gression up to v" — 39 with the present experimental 
accuracy will not reduce the uncertainties of Cq and Do 
significantly. This could be reached, however, even with 
the present distribution of data, provided the accuracy 
can be improved. Therefore, we believe that only high 
resolution Doppler free spectroscopic techniques should 
be used for further improvement on the Ca2 ground state 
PEC, especially at long internuclear distances. 
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FIG. 7: Distribution of the s-wave scattering length (a) and the position of the last bound level with respect to the asymptote 
(b) for the long range potentials obtained from the Monte Carlo simulation. 
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